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ABSTRACT 

Context. IGR J18410-0535 is a supergiant fast X-ray transients. This subclass of supergiant X-ray binaries typically undergoes few- 
hour-long outbursts reaching luminosities of 10^*- 10^^^ erg/s, the occurrence of which has been ascribed to the combined effect of the 
intense magnetic field and rotation of the compact object hosted in them and/or the presence of dense structures ("clumps") in the 
wind of their supergiant companion. 

Aims. IGR J18410-0535 was observed for 45 ks by XMM-Newton as part of a program designed to study the quiescent emission of 
supergiant fast X-ray transients and clarify the origin of their peculiar X-ray variability. 
Methods. We carried out an in-depth spectral and timing analysis of these XMM-Newton data. 

Results. IGR J18410-0535 underwent a bright X-ray flare that started about 5 ks after the beginning of the observation and lasted 
for ~ 15 ks. Thanks to the capabilities of the instruments on-board XMM-Newton, the whole event could be followed in great detail. 
The results of our analysis provide strong convincing evidence that the flare was produced by the accretion of matter from a massive 
clump onto the compact object hosted in this system. 

Conclusions. By assuming that the clump is spherical and moves at the same velocity as the homogeneous stellar wind, we estimate 
a mass and radius of Md==1.4xl0-^ g and i?d^8xlO" cm. These are in qualitative agreement with values expected from theoretical 
calculations. We found no evidence of pulsations at ~4.7 s after investigating coherent modulations in the range 3.5 ms-100 s. A 
reanalysis of the archival ASCA and Swift data of IGR J18410-0535, for which these pulsations were previously detected, revealed 
that they were likely to be due to a statistical fluctuation and an instrumental effect, respectively. 

Key words. X-rays: binaries - stars: individual IGR J18410-0535 - stars: neutron - X-rays: stars 
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1. Introduction than 1 hr and separated by 0.6 days. The peak fluxes were 2.0 

V ^ ,^,r, ^^^^ , .,,x,„.. „ ■ , ^ , and 9.5x10"" erg cm"^ s"' (2-10 keV), respectively. A time- 

^ IGR J18410-0535 (= AX J184L0-0536) is a member of the su- ^^^^j^^^ spectral analysis of the event revealed that the source 

O pergiant fast X-ray ttansients (SFXT), a subclass of supergiant ^-ray spectrum could be accurately modeUed with an absorbed 

^ X-ray binaries (sgHMXBs) that has atti-acted much attention in ^^^,.1^^ photon index r=l-2 and an absorption column 

the past few years because of their peculiar behavior in the X-ray of A^h=(3.2-7.2)x 10^^ cm-^. An iron line with a centroid 

domain (see, e.g. Walter & Zurita Heras 2007, for an updated list ^ 4 equivalent width (EW) of -0.2 keV 

of sources m this class . In contrast to previously known sgH- ^j^^ -^^^ ^ ^^^^ q^^- ^^-^^^^^ 

N^s which are nearly persistent in X-rays (see, e.g. |ChatiJ ^^^^^ pulsations were reported at a period of -4.7 s in only the 

120101 for a review), SFXTs usually undergo few-hour long out- ^ 9.4.9 keV energy band. The lowest X-ray flux of the source 

bursts reaching luminosities of lO^^'-lO^^ erg/s and spend most ^^^^^^^^ wifli A5CA was 2.0x10-'^ erg cm-^ s"' (2-10 keV). 
of their time life in quiescence , with typical luminosities of 
10^2-10^^ erg/s (Sidoh 2010; B ozzo et al.li2010l) . The presence 

of black-hole accretors in these sources cannot be completely IGR J1841( )-0535 was observed with Chandra for 20 ks on 

ruled out. However, the properties of their X-ray spectra in out- 2004 May 12 (|Halpern et alj|2004- This observation captured 

burst and quiescence and the detection of pulsations in the X-ray the source with an average flux of 4.2x10 erg cm ^ s ' (0.5- 

emission of some SFXTs, led to the conclusion that (at least) 10 keV), and the corresponding X-ray spectrum could be de- 

most of fliem shou ld host neutron-star (NS) accretors (see e.g., scribed well by an absorbed power-law model with r= 1 .35+0.30 

iBozzo et al ."2008b, and references therein). and A^h=(6.1 + 1.0)x10^^ cm ^. This observation also provided 

IGRJ18410-0535 was discovered with ASCA in 1994 an improved position of the source at aj2ooo=18HF"0:54 and 

jBamba et al.1 120011) . During the discovery observation, the 5j2000=-05°35'46'.'8 (nominal Chandra position accuracy 0.6'.'), 

source displayed two bright flares with rising times shorter and permitted the identification of the optical and infrared coun- 
terparts of the source a s a B 1 lb supergian t star at an estimated 



Send offprint requests to: E. Bozzo distance of 3.2^j 9 kpc (iNespoli et al.ll2008l) . 
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The source was also captured undergoing several few- 
hour-long outbursts with MAXI (Negoro et al. 201 0) and 
INTEGRAL ( 1-7 hr; [Rod riguez et al. 2004; Sguer aet all 
120061: IWalter & Zurita HerasI l2007h . The highest flux mea- 
sured by IN TEGRAL du r ing th ese events was 120 mCrab 
(20-80 keV; ISguera et"al] 120061) . corresponding to roughly 
1.8x10"^ erg cm"^ s"'. The source high energy spectrum (20-80 
ke V) while in outburst could be described using a hot blackbody 
model (BB) with kT-%-9 keV. Outside the outbursts, the upper 
limit (Icr c.l.) to the average hard X-ray flux of the source was 
estimated to be 1 mCra b (18-60 keV, corre sponding to roughly 
1.2x10 " erg cm 2 s''. lFilippova et al.ll2004l) . 

IGRJ18410-0535 was detected a few times in out- 
burst by 5w//f/BAT, and only once Swift performed a slew 
on the sou rce to point it with t he narrow field inst ru- 
ment XRT dde Pasquale et all 1201 Ot iRomano et al. I l20T0allbh . 
On that occasion, the 5w//fXRT spectrum was closely de- 
scribed by an absorbed power-law model with r=0.7;^Q^ and 
A^H=3;^jXlO^^ cm"^. The average 2-10 keV unabsorbed flux 
was 7x10""^ erg cm"^ s"'. No evidence for pulsations was 
found. IGRJ18410-0535 was not detected in a high X-ray 
activity phase during the one-year monitoring (total net on- 
source exposu r e time 96.5 ks) with Swift jXKl carried out by 
iRomano et alJ (l2009h . In the data acquired during this moni- 
toring, the authors identified four different states for the X-ray 
emission of IGR J18410-0535: the "high", "medium", "low", 
and "very low" states. In the high state, the source X-ray flux was 
8.0x10"" erg cm"^ s"' (2-10 keV) and the corresponding spec- 
trum is closely described by an absorbed power-law model with 
r=l. 1+0.1 and A^H=(2.5±0.3)xl022 cm^^ j^e medium state 
was characterized by a flux of 3.4x 10"" erg cm"^ s"' , a power- 
law photon index of 1.3+0.2, and A^h=(3.5+0.5)x10^^ cm"^. 
For the low and very low states, the authors estimated fluxes 
of 1.1x10"" erg cm"^ s"' and 6.0x10"'^ erg cm"^ s"', re- 
spectively. The corresponding power-law photon indices and 
absorption column densities were r= 1.5+0.1, 0.6+ 0.4 and 
A?W3.5±0.5)xl022 cm (0.6+0.4)xl022 cm l Sid oli et alj 
(|2008|) also reported the detection of pulsations at ~4.7 s from 
IGRJ18410-0535 in some of the Swift IXm data. A possible 
association between IGR J18410-0535 and the tr ansient MeV 
EGRE T source 3EG J1837-0423 was suggested bv lSguera et alj 
(I2009h . 

In this paper, we report on a 45 ks long observation of the 
source with XMM-Newton. This observation was obtained as 
part of our progra m aimed at s tudying the quiescent emission 
of SFXTs (see also lBozzo et "an i2010. 2009. 2008b). During the 
XMM-Newton observation, IGRJ18410-0535 was caught under- 
going a bright X-ray flare, lasting roughly 15 ks. The source 
reached a maximum flux of 3.4xl0""^erg cm"^ s"^ (1-10 keV) 
at the peak of the flare and then decayed to a very low quiescent 
level (8.8xlO"'"^erg cm"^ s"'). The total dynamic range in the X- 
ray flux was thus >4xl0^. In Sect. |2] we give the details of our 
data analysis and provide all the results in Sect. 12.21 and l2. 31 In 
the XMM-Newton data, no pulsation at 4.7 s could be detected, 
and the derived upper limits were significantly tighter than the 
pulsation amplitude reported previously with ASCA and Swift. 
Motivated by these results, we also reanalyzed all the ASCA and 
Swift data of the source and show that the previously reported 
4.7 s pulsations were probably due a statistical fluctuation {ASCA 
) and an instrumental effect (Swift ), respectively (Sects. [3]and|4|. 
Our discussion and conclusions are reported in Sects. |5]and|6] 
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Fig. 1. XMM-Newton Epic-pn lightcurve of the observation of 
IGR J18410-0535 (0.3-12 keV, not corrected for background and 
pile-up). The start time is MJD 55270.5896 (bin time 30 s). 
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Fig. 2. Rise of the flare observed during the XMM-Newton ob- 
servation of IGR J18410-0535 with the Epic-pn camera. The 
start time of the lightcurve is the same as that in Fig. [1] and 
the binning time is 150 s. The upper panel shows the source 
lightcurve in the soft energy band (0.3-4 keV), the middle panel 
in the hard energy band (4-12 keV), and the bottom panel the 
hardness ratio (HR, defined as the ratio of the source count-rate 
in the soft to hard energy band versus time.) 



2. XMM-Newton data analysis 

XMM-Newton (iJansen etal.ll200lh observed IGR J18410-0535 
from 55270.5484 MJD to 55271.0 923 MJD (e xposure time 
~45 ks). The Epic-pn camera (Striider " ;t alJl200ll) was operated 
in full frame , while the Epic-MOSl and Epic-MOS2 cameras 
([Turner et al.ll2001i) operated in small window and timing mode, 
respectively. This particular set-up of the instruments was cho- 
sen in order to have the capabilities required to study the source 
over a wide range of X-ray flux. We processed the XMM-Newton 
observation data files with the two pipelines epproc and emproc 
in order to produce Epic-pn and Epic-MOS cleaned event files, 
respectively (SAS v.10.0.1). No time intervals were found to be 
affected by a high background and we thus retained for the sub- 
sequent analysis the entire exposure time available for the three 
instruments. We used the energy range 0. 15-15 keV for the Epic- 
pn and 0.3-12 keV for the two Epic-MOS cameras. IGR J18410- 
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Fig. 3. Details of the Epic-pn lightcurve during the drop in count 
rate that occurred around f=8500-9500 s in Fig. [1] (f=0 would 
correspond to the the start time of the lightcurve in Fig.[Tl time 
bin is 30 s). The upper (lower) panel shows the source lightcurve 
in the 0.3-4 keV (4-12 keV) energy band. 
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Fig. 4. Final part of the decay after the flare observed by XMM- 
Newton. The upper (lower) panel shows the source lightcurve in 
the 0.3-4 keV (4-12 keV) energy band. We marked with vertical 
dashed lines the time intervals during which major changes in 
the source spectrum occurred (spectra K,L,M,N,0 in Table[T]and 
Fig.EJ. 



0535 displayed a large dynamical range in the X-ray flux during 
the XMM-Newton observation (see Fig.[TJ; therefore the source 
and background extraction regions were chosen in the different 
time intervals so as to maximize the signal-to-noise ration (S/N) 
of the data (see next section). During the periods of greatest X- 
ray emission from the source, the Epic-pn data suffered signif- 
icant pile-up. We accounted for this problem by using in these 
cases annular extraction regions for the source in which the ap- 
propriate innermost portion of the instrument point spread func- 
tion (PSF) was removecQ. We checked that the removal of pile- 
up was effective by using the SAS tool epatplot and comparing 
the results from the Epic-pn camera with those obtained with the 
Epic-MOS 1 and Epic-MOS2. The effect of pile-up in these cases 
was much reduced since the cameras were operated in small 
window and timing mode, respectively. Where required, we cor- 



' See also http://xn im.esac.esa.int/extemal/xmm_user_supportl /docu- 
mentation/uhb/index.html. 
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Fig. 5. The same as Fig. [T] but here we used an adaptive rebin- 
ning with S/N=5. The upper panel shows the source lightcurve 
in the 0.3-4 keV energy band, the middle panel the source 
lightcurve in the 4-12 keV energy band, and the lower panel the 
hardness ratio versus time (defined as in Fig.|2|. 

rected Epic images to remove the out-of-time (OoT) event^ and 
checked that this problem did not significantly affect the spectra. 
All the Epic lightcurves were barycentered by using the SAS 
tool BARYCORR, and then corrected for instrumental vignetting, 
dead time, and PSF losses by using the tool epiclccorr. Given 
the lower count rate and S/N of the Epic-MOS 1 and Epic-MOS2 
cameras compared to those obtained from the Epic-pn (espe- 
cially outside the flare), the former instruments did not con- 
tribute significantly to the spectral analysis. Therefore, we report 
in the following sections and in Fig. [1] only the results obtained 
from the Epic-pn camera. 

2.1. Lightcurve analysis 

During the XMM-Newton observation, IGR J18410-0535 under- 
went a bright X-ray flare that started about 5 ks after the begin- 
ning of the observation and lasted for ~15 ks. The lightcurve of 
the entire XMM-Newton observation in the 0.3-12 keV energy 
band is shown in Fig.[T] 

The lowest source count-rate was recorded during the latest 
20 ks of observation at a level of (4.0+1.9)xl0"-' cts/s (estimated 
from a 5000 s binned lightcurve that had background subtracted; 
all uncertainties in this paper are given at 90% c.l., unless other- 
wise indicated). During the peak of the flare, the corresponding 
highest count-rate was 31.2±1.0 cts/s (determined from a 30 s 
binned lightcurve, which was background subtracted and cor- 
rected for pile-up). In Fig. |2] we show an enlargement of the 
lightcurve during the first 5 ks of the XMM-Newton observation. 
Here, a sudden rise of the source X-ray emission was recorded 
by the Epic-pn camera ~3 ks after the beginning of the obser- 
vation. In the following ~2 ks, the source reached the highest 
count-rate, and then began to decrease at f~6000 s. Interestingly 
a drop in the source count rate was recorded around r=8500 s in 
both the soft and hard energy band. A zoom in of this part of the 
lightcurve is provided in Fig. [3] A detailed lightcurve of the final 
part of the flare decay is shown in Fig.|4] We also show in Fig.|5] 
the source lightcurves in the 0.3-4 keV and 4-12 keV energy 
bands, where an adaptive rebinning of the data was used in order 

^ See |http://xmin.esa.int/sas/current/documentation/threads/| 
EPIC_OoT.shtml. 
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IGRJ18410-0535 (Epic-pn) 
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Fig. 6. Hardness-intensity diagram of IGR J18410-0535 ob- 
tained by using the XMM-Newton Epic-pn observation. The 
hardness ratio of the source is defined as in Fig.|2l 




Fig. 7. Spectral parameters of IGR J18410-0535 measured dur- 
ing the different time intervals reported in Table [T] All the errors 
are at 90% c.l. Arrows indicate 90% c.l. upper limits. 



to achieve in each time bin S/N=5 and estimate the correspond- 
ing source hardness ratio (HR). The same lightcurves and adap- 
tive rebinning were used to produce the hardn ess ration ver sus 
intensity diagram shown in Fig|6](see also Boz zo et alj|2010l) . 

From Figs. |5] and|6] we can see that the HR of the source 
underwent dramatic changes in time; a non-monotonic trend ap- 
pears when the HR is plotted as a function of the total source 
intensity. Given these findings, we performed in Sect. l2.2l below 
a time-resolved spectral analysis of the data. 

2.2. Spectral analysis 

To search for spectral changes, we first accumulated X-ray spec- 
tra on time intervals of ~few hundreds to ~few thousands sec- 
onds (depending on the source count-rate), following the trend of 
the HR (see Fig. |5]l. Nearby time intervals were then combined 
in order to extract a lower number of spectra that could max- 
imize the evidence for changes in the properties of the source 
X-ray emission. Our optimal separation of the time intervals is 
reported in Table [T] All these spectra were rebinned to have at 
least 15, 20, or 25 photons per bin (depending on the source in- 
tensity) and prevent an oversampling of the energy resolution 
of the instruments by more than a factor of three. Spectra with 
lower quality statistics were rebinned t o ha ve at least 5 photons 
per bin and were fit using C-statistics (lCashlll97 9). For the spec- 
tral fits, we used an absorbed power-law modefl None of the 
spectra could be convincingly better fit with a blackbody (BB) 
or a cut-ofF power-law model (cutoffpl in Xspec). In particular, 
the latter provided in most cases peculiar values of the model 
parameters (negative power-law photon indices and values of 
the high energy cut-ofF largely inconsistent between the differ- 
ent spectra). A plot of the evolution of the spectral parameters 
of the source in the XMM-Newton observation is provided in 
Fig. I2I The power-law photon index remained fairly constant 
up to 19000 s, whereas the absorption column density under- 
went dramatic changes (Fig. |7]). In particular, we measured a 
large increase (factor of ~10) in the absorption column density 
from the beginning of the observation up to the peak of the flare 
(A^H~2xlO^^ cm"^ at f^5000 s). A rapid drop in the A^h then 

^ We used the model phabs*pow in xspec. 
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Fig. 8. XMM-Newton Epic-pn spectra of IGR J18410-0535 ex- 
tracted during the latest part of the flare decay (see Table [1] and 
Fig. nil. All the spectra were fit with a simple absorbed power- 
law model (the residuals from these fits are shown for each of 
the five spectra). 



occurred to ~10^^ cm"^. At the peak of the flare, we also de- 
tected a significant iron line with an energy of ~6.6 keV and 
an EW~60 eV (see Table[T]). The absorption column density re- 
mained virtually constant at ~10^^ cm"^ for about 4 ks after the 
source reached the highest X-ray emission level; a rapid rise in 
column density then occuiTed until f^l9000 s. To test the signif- 
icance of the measured variation in the absorption column den- 
sity, we report in Fig.|9]the A^h-F parameters confidence contours 
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Table 1. Best-fit parameters of IGR J18410-0535 during different time intervals of the XMM-Newton observation (see Fig. [1] and 
the note below the table). 



Interval 


Tstart 


Tstop 


Exp. 


Wh 


r 


^line 








^unabs 








(s) 


(s) 


(ks) 






(keV) 


(keV) 


(erg cm" 


2s-l) 


(erg cm- 


2 s-') 


(C-statistics/d.o.f.) 


A 





3200 


2.7 


2 6+2-5 


1 3+''' 


— 


— 


3.2x10- 


13 


4x10-'' 




(12.4/11) 


B 


3200 


4670 


1.2 


8.0+5-« 


1.1±0.9 


— 


— 


1.7x10- 


12 


2.8x10- 


12 


(14.4/18) 


C 


4670 


5020 


0.35 


19.3;|3 


2.2±0.7 


— 


— 


1.0x10- 


11 


5.0x10- 


11 


(51.0/47) 


D 


5020 


5220 


0.17 


19.8!80 


1.9±0.7 


— 


— 


6.0x10- 


11 


2.4x10- 


10 


0.8/12 


E 


5220 


6320 


0.96 


10.8!«;^ 


1.1±0.1 


6.56±0.05 


0.06 


3.2x10- 


10 


6.0x10- 


10 


1.0/149 


F 


6320 


6820 


0.44 


10.9!|;2 


1.1±0.1 


— 


— 


2.9x10- 


10 


5.3x10- 


10 


0.9/109 


G 


6820 


7620 


0.7 




1.4±0.2 


— 


— 


2.2x10- 


10 


4.4x10- 


10 


1.1/108 


H 


7620 


8520 


0.8 




1.5±0.2 


— 


— 


1.3x10- 


10 


2.9x10- 


10 


1.1/73 


I 


8520 


9520 


0.9 


15 5+'-'' 


1.6±0.2 


6.32±0.05 


0.10 


6.2x10- 


11 


1.8x10- 


10 


0.8/107 


J 


9520 


10420 


0.8 


14 3+' -'* 


1.2±0.2 


— 


— 


8.4x10- 


11 


1.9x10- 


10 


0.8/110 


K 


10420 


11520 


1.0 


17.0±1.5 


1.3±0.2 






6.5x10- 


1 1 


1.6x10- 


10 


1.1/114 


L 


11520 


14520 


2.6 


23.3±2.0 


1.6±0.2 






2.3x10- 


11 


8.5x10- 


11 


1.0/102 


M 


14520 


16620 


1.8 


28.0!5-0 


1.7±0.4 






8.2x10- 


12 


3.3x10- 


11 


1.2/58 


N 


16620 


18870 


1.8 


4T 4+10.7 


2.0±0.6 






5.1x10- 


12 


5.3x10- 


11 


1.4/36 





18870 


42840 


20.8 


<0.5 


-0.9±0.3 


6.36±0.06 


1.1 


2.2x10- 


13 


2.2x10- 


13 


1.0/13 


Oi 


18870 


22820 


3.4 


<8.0 


■l-**-0.4 


6.33±0.06 


1.2 


9.7x10- 


13 


9.8x10- 


13 


0.8/5 


02 


22820 


42840 


17.4 


<1.4 


3+''-5 






8.8x10- 


14 


9.1x10- 


14 


(12.5/10) 



Notes. The continuum spectral model is an absorbed power law. Here, A'h is the absorption column density, Y is the power-law photon index, Snne 
is the energy of the centroid of the iron line, and £Wii„c the corresponding equivalent width. Fobs (Funabs) is the absorbed (unabsorbed) flux in the 
1-10 keV band. In the last column we report the value of the;^'2^j/d.o.f or the corresponding value of the C-statistics/d.o.f, for each spectrum. 

for some relevant spectra. In the upper panel of this figure, we 
show the parameter contours obtained by summing the spectra 
from intervals C and D and A and B, which displayed very simi- 
lar values of the power-law photon index and absorption column 
density before and during the rise of the flare, respectively (see 
Table|T]i. The contours correspond to 68%, 90%, and 99% c.l.. A 
very significant (>3 cr) increase in the absorption column den- 
sity at the onset of the flare is clearly seen. The lower panel of 
Fig. |9] shows the case of the J and N spectra: the increase in 
the absorption column density toward the end of the flare is also 
highly significant. During the apparent drop in count rate around 
f=8500-9500 s, no particular change in the spectral continuum 
is seen, but a significant iron line appears with a centroid energy 
of ~6.4 keV (see spectrum I in Table [TJ throughout this paper, 
we use narrow lines for the spectral fits and fixed their width to 
in XsPEc). 

The properties of the X-ray emission of the source changed 
sharply around 19000 s: at this time, the source underwent 
a further decrease in the X-ray flux, and simultaneously its 
spectrum flattened and a prominent iron emission line appeared 
around 6.4 keV. We show in Fig. [8] the detail of the evolution 
of the source spectrum with time during the latest 30 ks of the 
XMM-Newton observation. A fit to the spectrum from the O in- 
terval with a simple absorbed power-law model gave unaccept- 
able results (x^^Jd.o.f.=2.\l\5). Adding an iron line with an en- 
ergy centroid of £'iine=6.36;'^Qp5 keV significantly improved the 
fit ix^^Jd.o.f.-l.Qliy). The normalization and equivalent width 
of the line were (2.7+1.0)xl0"* and 1.1 keV, respectively. From 
this spectrum, only an upper limit on A^h could be obtained (see 
Table [T]|; this was significantly lower than that expected in the 
direct ion of the source (~ 1.8x10^^ cm'^, iDickev & LockmanI 
119901) . Even though this model gave a statistically acceptable 



fit to the spectrum O, the very low value of the absorption col- 
umn density and the unphysical negative power-law photon in- 
dex suggest that other spectral models cannot be excluded. As 
we discuss in Sect. |5j a more physical interpretation of the X- 
ray emission of the source in this time interval is based on 
a modefl comprising two power-law components affected by 
two different absorption column densities plus the iron line, 
i.e. PHABsl*(powl-i-GAUss+PHABs2*pow2) in XsPEC. In the fit, we 
fixed phabsl to the Galactic value of the absorption column den- 
sity and constrained the photon indices of the two power-laws 
to be the same. The fit gave ;i'2^yd.o.f. =0.8/12, r=0.9+';^ (in 
agreement with that measured during the other time intervals in 
Table[B, £iine=6.36+0.06 keV, £W=0.7 keV, and a value for the 
second absorption component of (85;'^33)xlO^^ cm"^. In Fig.fTOl 
we report the unfolded spectrum of interval O together with the 
different spectral component used in this model. 

In Fig. [TT] we also show the contours plots of the iron line 
centroid energy versus its normalization for the three detections 
in the E, I, and O spectra. We find marginal evidence for a de- 
crease in the line centroid with time, from the peak of the flare 
(spectrum E) to the quiescent level (spectrum O). No other con- 
vincing evidence of the iron line could be found in the spectra 
extracted from different time intervals or from merging them in 
different combinations. The interpretation of all the above find- 
ings is discussed in Sect. 15. II 

We also performed a separate spectral analysis of the source 
X-ray emission during the first ~3 ks and the last ~17 ks of the 
time interval O in order to investigate the cause of the drop in the 

We also performed a fit using the reflection model pexrav in Xspec, 
but the poor statistics of the O spectrum did not permit to derive satis- 
fiable constraints on the model parameters. We thus do not discuss this 
model in further details. 
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IGRJ 18410-0535 Contour plots (Epic-pn) 




Fig. 9. Upper panel: Contour plots for the spectral parameters 
measured from the time intervals A+B and C+D (see Table [1]). 
The contours correspond to the 68%, 90%, and 99% c.l. Lower 
panel: Same as above but for the J and N spectra in Table [T] 

IGRJI84I0-0535 (Epic-pn, Unfolded sptjctriitn "O") 



> 




Fig. 10. Unfolded Epic-pn spectrum of IGR J18410-0535 during 
the interval O in Table [1] The best-fit model in Xspec here is 
phabsl*(powl-i-gauss+phabs2*pow2). The value of phabsl was 
fixed in the fit to 1.8x10^^ cm"^ and the photon index of the 
two power laws was forced to be the same (see Sec. l2.2l for de- 
tails). The dashed line corresponds to the less absorbed power- 
law component in the fit, and the dot-dashed line to the more 
absorbed one. The dotted line indicates the iron line component 
at ~6.4 keV. The bottom panel shows the residuals from the fit. 



source count-rate visible in Fig. |4] around f =22820. We indicate 
these two additional spectra with Oi (from f=18870 to f =22820) 
and O2 (from f=22820 to f=42840) in Table[I] 

The Oi spectrum comprised 151 counts and was dominated 
by the presence of the iron line. A fit with a simple absorbed 
power-law did not provide an acceptable fit (x^^Jd.Q.f.-2.9ll), 
whereas adding a Gaussian line to the spectral model sig- 
nificantly improved the fit (x^^Jd.o.f.-Q.%l5). We measured a 
power-law photon index of r=-1.8^y g, and a centroid energy 
for the line of 6.33^^ jj^ keV. The normalization of the line was 
(1.1 +0.4)x 10"-'' and the estimated EW 1 .2 keV (compatible with 
that found by using the entire time interval O). Only an upper 
limit to the absorption column density of A^h<8.0x10^^ cm"^ 
could be obtained. 

Given the measured negative value of the parameter F, we 
also fit the spectrum Oi by using the double power-law model 
discussed above and checked the consistency of the different 
model parameters with those determined before for the spectrum 
O. This gave F=0.2^[ ° and an absorption column density for the 
most extinguished power-law component of 82;^4^xl0^^ cm"^ 
(the absorption column density of the other component was 
fixed to the Galactic value). The centroid energy of the line 
and its EW was found to agree with those reported before 
for the spectrum O. The estimated 1-10 keV X-ray flux was 
9.7x10"''' erg cm"^ s"'. 

During the time interval O2, XMM-Newton recorded only 
83 counts from the source. The corresponding spectrum could 
be closely fit with a simple absorbed power-law model (see 
Table [U. The quality of the fit did not change significantly (C- 
statistics 14.5/11) when the power-law photon index was con- 
strained to be equal to that measured for the O spectrum (F= 1 .0); 
the corresponding value of the absorption column density was 
A?h=(1-6^[ q)x10^^ cm"^, compatible with the expected Galactic 
value in the direction of the source. Given these results, we did 
not attempt to fit this spectrum with a double power-law model. 

Adding to the spectral model an iron line with an energy 
fixed at £'iine=6.36 keV provided an upper limit (90% c.l.) to 
its normalization and EW of 1.8x10"^ and 1.9 keV, respectively. 
While the upper limit to the EW is consistent with that estimated 
for the O spectrum, the normalization is slightly lower However, 
it is still compatible with the value expected because of the de- 
crease in the flux in the two spectra. 

All these results suggest that the drop in the count-rate no- 
ticed around f~23000 s was accompanied by the occurrence of a 
further spectral change in the X-ray emission from the source. 
The spectrum Oi could not be described well by assuming a 
simple absorbed power-law model and the addition of a second 
more absorbed component (plus the iron line) provided a more 
reasonable fit to the data. In contrast, the spectrum O2 could be 
well described by using a single relatively low absorbed power- 
law component. We show the two spectra Oi and O2 in Fig. [12] 
and discuss their interpretation further in Sect. |5] 

We note finally that the relatively low quality statistics of 
spectra A and O2 and large errors aff'ecting the values of their 
best-fit parameters prevent us from making a clear statement 
about any change in the spectral properties of the source before 
and after the occurrence of the flare. The XMM-Newton data re- 
vealed, however, that the flux of the source before the onset of 
the event (time interval A) was a factor ~3.5 higher that that es- 
timated during the last part of the observation (time interval O2). 
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IGRJ1841()-()535 (Epic-pn spectrum E) IGRJ1841.0-0535 (Epic-pn spectrum I) IGRJ18410-0535 (Epic-pn spectrum O) 




5^-5 2 5 

Energy (keV) Energy (keV) Energy (keV) 



Fig. 11. Top panels: Contour plots of the iron line centroid energy vs. normalization measured from the spectrum E (left), I (middle), 
and O (right). The contours correspond to the 68%, 90%, and 99% c.l. There is a marginal indication of a change in the centroid 
energy of the line between the time interval E and the other two. Bottom panels: The ratio of the data to the model of the spectra 
extracted during the time interval E (left), I (middle), and O (right). All these spectra were fit by using only an absorbed power-law 
model in order to highlight the presence of the iron line in the residuals from the fit. 



2.3. Timing analysis 



IGRJ 184 10-0535 (Epic-pn, spectra O, and O,) 



IGRJ18410-0535 was report ed in previous studies to emit pul- 
sation s at a period of ~4.7 s dBamba et alj|2001l : iRomano et al.l 
l2009l see also Sect. [TJ. We carried out an in-depth search 
for coherent modulations in data collected with both the Epic- 
pn and Epic-MOS cameras by extracting barycentered source 
and background event lists (we used the most accurately de- 
termined source position; see Sect. [T] and I2.3l l. We applied to 
the lists of barycentered photon arrival times in the 0.3-12 keV 
energ y ra nge the power-spectrum search algorithm developed 
by ITsrael & Stellal (Il996h . This method is optimized to search 
for periodicities in "coloured" power spectrum components and 
derive upper limits if no signal is detected. Pulsations were 
searched for by using both the photon arrival times of the en- 
tire XMM-Newton observation and in the different time intervals 
selected for the spectral analysis. We did not find any significant 
indication of pulsations. As an example. Fig. \13\ shows the re- 
sults obtained by applying the method of llsrae l & Stella ([19961) 
to the list of photon arrival times extracted from the entire XMM- 
Newton observation. In Tablej2l we also report all the 3 cr upper 
limits on the pulsed fraction^ and pulsation frequencies we in- 
vestigated. We show in this table the upper limit estimated by 
using XMM-Newton data accumulated during the entire observa- 
tion (A+B+C+D+E+F-hG+H+I+J+K+L+M+N+O), the flare 
event (including rise and decay, C-i-D+E+F+G+H), the rise of 
the flare (C-hD), the top of the flare (E), the beginning of the flare 
decay (F-i-G-i-H), the last part of the flare decay (J-hK-i-L), and 
the time interval in which the higher absorption was measured 
(M-i-N). This selection of intervals was carried out to investigate 
whether pulsations might have been present only during part of 
the event observed by XMM-Newton. 

For all the time intervals reported in Table |2l we also per- 
formed additional searches of pulsations by separating the lists 
of photon time arrival times in the soft (0.3-4 keV) and hard 




^ Here we define the pulsed fraction as the semi-amplitude of the si- 
nusoid divided by the source average count rate, (/max-/min)/(Anax+/inin)- 



Energy (keV) 

Fig. 12. The two spectra Oi and O2 extracted during the fist 
~3 ks and the last ~17 ks of the time interval O in Table [1] To 
facilitate the visual comparison, we used in this figure a group- 
ing of 5 photons per bin for both spectra. We show the fit ob- 
tained with a simple power-law model to highlight the iron line 
at -6.4 keV. 



(4-12 keV) energy bands and by combining strictly simultane- 
ous EPIC-pn and MOS event lists (in this case, a common time 
resolution of 0.3-s was adopted). No significant improvement to 
the upper limits indicated in Table |2] could be obtained. Even 
though for a number of time intervals considered in Table |2] the 
derived 3 o" upper limit to the pulsed fraction was relatively high 
(~30-50%), the entire exposure time of the XMM-Newton ob- 
servation provided relatively tight constraints on the presence of 
pulsations from IGR J18410-0535 compared to those reported 
previously (see Sect.[T]i. Motivated by these findings, we reana- 
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Table 2. Three cr c.l. upper limits on the pulsed fraction of IGR J18410-0535 (see notes below the table). 



XMM^Newton (0.3-12 keV) 


Time interval 






Instrument (time resolution) 






r\\ \ 1 j.^ rns } 




ivnj^ 1 yyj.j s ) 


iVlUoZ i . / J mS ) 


A+B+C+D+E+F+G+H+I+J+K+L+M+N+O 


6-10% (0.02-6.8 


Hz) 


10-20% (0.01-1.6 Hz) 


15-30% (0.02-286 Hz) 


C+D+E+F+G+H 


10-30% (0.04-6. 


8 Hz) 


10-30% (0.01-1.6 Hz) 


10-30% (0.03-286 Hz) 


C+D 


30-50% (0.05-6. 


8 Hz) 


50-70% (0.1-1.6 Hz) 


60-80% (0.06-286 Hz) 


E 


20-35% (0.05-6. 


8 Hz) 


20-35% (0.03-1.6 Hz) 


20-70% (0.1-286 Hz) 


F+G+H 


15-20% (0.08-6. 


8 Hz) 


20-30% (0.02-1.6 Hz) 


20-30% (0.05-286 Hz) 


J+K+L 


10-20% (0.02-6. 


8 Hz) 


20-40% (0.01-1.6 Hz) 


35-60% (0.02-286 Hz) 


M+N 


40-60% (0.02-6. 


8 Hz) 


75-100% (0.02-1.6 Hz) 


80-100% (0.02-286 Hz) 


ASCA (0.7-10 keV) 


GIS (0.5 s) 


Obs.l+Obs.2 


30-45% (0.01-1 


Hz) 







Notes. The three a upper limits on the pulsed fraction of IGR Jl 8410-0535 are determined for the three Epic cameras in the different time intervals 
reported in Table[T] The results of a similar analysis carried out on the ASCA data are also reported (see Sect. l3.1l for the discussion of the ASCA 
data analysis). 




10-* 10-3 01 1 

Frequency 



Fig. 13. Power spectrum produced using data from the entire 
XMM-Newton observations (0.3-12 keV). The upper continuous 
line represents the power threshold for detection of periodicity at 
a 3 cr C.L, according to the method described bv llsrael & Stellal 
We also show in the bottom part of the plot a second con- 
tinuous curve that represents the upper limit calculated with the 
same method, on the presence of pulsations from IGR J 184 10- 
0535 as a function of the frequency. The two dashed lines repre- 
sent the 5% and 10% upper limit levels on that curve. The most 
stringent upper limit we could provide with this method is 6-10% 
for frequencies in the range 0.02-6.8 Hz. 



lyzed all the published ASCA and Swift data for which detections 
of pulsations at ~4.7 s were reported. 



Table ^.ASCA observations log of IGR JI84I0-0535. 



Instmment 


Obsei'vation 


Start time 


Time elapsed 


Expo .sure 


ASCA(l) 


61009100 


1994-12-04X15:48:13 


4312.0 


2112.0 


ASCA (2) 


57026000 


1999-10-03105:47:32 


65484.0 


37373.6 



Imaging Spectrometers ("SIS. lBurke et a l."l993'). but included in 
that of Gas Imaging Spectrometers (GIS, Ohashi et al. 199^. 
The GISs were operated in the nominal PH mode with a time 
resolution of 62.5 ms for data taken in high bit-rate and 0.5 s for 
those taken in medium bit-rate. For our analysis, we used only 
GIS data obtained when the satellite was either outside the South 
Atlantic Anomaly and low cut-off rigidity regions (> 6 GV), or 
when the target elevation angle was >5°. Particle events were re- 
moved b ased on ris e-time discrimination with the gisclean task 
(Ohashi e t al.lll99^ . After this screening, the total available ex- 
posure times of Obs. 1 and 2 were 2 ks and 36 ks, respectively. To 
improve the quality of the statistics, we combined data from the 
GIS-2 and GIS-3 detectors. ASCA events were extracted from a 
circular region of 3 arcmin radius centered on the best known 
source position (Halpern et al. 2004). The background was ex- 
tracted from an annular region around the source position with 
inner radius of 4.3 arcmin and outer radius of 6.8 arcmin. This 
annular region is free from field point sources and is unaffected 
by the extended diffuse emission from the nearby supernova 
remnant G26. 6 -0.1 (Bamba et al. 2003). The inner radius of the 
background extraction radius was chosen so as to avoid contami- 
nation due to the tail of the ASCA PSF. Photon aiTival times were 
all converted to the Solar System barycenter (SSB ) with the task 
TiMECONV using the Chandra position reported by 'Halpe rn et al.l 
(2004). The log of the ASCA observations of IGR J18410-0535 
is given in Table [3] 



3. >tSCiA data analysis 

IGR J18410-0535 was observed by ASCA for the first time on 
1994 April 12 (MJD 49454.677-49454.701), during the survey 
of the Scutum arm region (hereafter Obs.l). A follow-up obser- 
vation was carried out on 1999 October 3-4 (MJD 51454.252- 
51455.463, hereafter Obs.2). On both occasions, IGRJ18410- 
0535 was outside the field of view (FOV) of the Solid-state 



3.1. Timing analysis 

In Fig. [141 we show the background-subtracted light curve (top 
panel) of IGR J18410-0535 in the 0. 7 - 10 keV e n ergy band. 
According to the results reported by Bamba et al.' (2001), we 
note that the source was characterized by a variable intensity, and 
that a 5 hr-long X-ray flare is visible at the end of Obs.2 (starting 
from MJD 51455.3). During this flare, the intensity of the source 
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Fig. 14. rop panel: ASCA light curve of IGR J18410-0535 in the 
0.7-10 keV energy band. Bottom panel: Hardness ratio of the 4- 
10 and 1-4 keV GIS light curves. The vertical dashed lines mark 
the flare episode starting at MJD 51455.3. The hashed region 
marks the separation between Obs.l and Obs.2. A and B denote 
the two data set selected for the timing analysis. 



o 

Q. 




Frequency (Hz) 

Fig. 15. ASCA power specti'um of IGR J18410-0535 obtained 
by using events in the 1.9-4.9 keV energy range and a time res- 
olution of 0.5 s. The arrow marks the apparent detection of pul- 
sations at 0.21 1 Hz. As we discuss in the text, this excess in the 
power-spectrum is not statistically significant if the correct num- 
ber of trials is taken into account (see Sect. 13. 11 1. 



increased by a factor of ~5. After -14 ks, IGR J18410-0535 re- 
turned to an X-ray emission level comparable to that observed 
before the flare. The bottom panel of Fig [141 shows the evolution 
of the hardness ratio computed by using the lightcurves extracted 
in the 0.7 -2 keV and 2-10 keV energy bands. In Obs.2, the HR 
showed a remarkable increase (from ~0.5 to ~2) with the source 
intensity. 

During the flare in Obs.2, iBamba et al.l (1200 lb reported on 
the detection of coherent pulsations by carrying out a fast Fourier 
transform (FFT) of the barycentered photons-event arrival times 
selected in the 1.9 - 4.9 keV energy range (energy channels 
161 - 415) and extracted from a 3' circular region centered on 
the best known source position. An excess in the power spec- 
trum (0.02-1 Hz) of -40 was visible at -0.211 Hz (their most 
accurately determined spin period was P=4. 7394+0. 0008 s, see 
also Sect. [T]i. Taking into account the total number of indepen- 
dent frequency trials, this excess corresponds to 4.4 (Gaussian) 
standard deviations. 

We revisited the data from Obs.2 performing a power spec- 
trum analysis. The source events were collected during the flare 
with a time bin of 0.5 s (the time resolution of medium bit- 
rate ASCA data) and extracted from a 3' radius circular region 
around the source. We investigated frequencies in the range 0.02- 
1 Hz with a resolution of 3.7x10"^ Hz, first by using the entire 
energy range of the GIS (i.e. 0.7 - 10.0 keV) and then only the 
events in the 1.9-4.9 keV and 4.9-10.0 keV energy intervals. 

In the first and the latter cases, the power spectrum did not 
reveal any significant deviations from a statistically flat distribu- 
tion. The corresponding estimated upper limit on the pulsed frac- 
tion of the source at different frequencies is reported in Table |2] 
(we used the method described in Sect. I2.3l in order to facilitate 
the comparison with the XMM-Newton results). The power spec- 
trum realized with the events in the energy interval 1.9-4.9 keV 
showed a peak at ~0.21 1 Hz with a power of ~33 (see FigfTSl). If 
all the inspected independent frequencies and the ad-hoc selec- 
tion in energy are taken into account to estimate the total effec- 
tive number of trials, this peak corresponds to a detection signif- 
icance of <3 (Gaussian) standard deviations. 



A signicance of ~4.4 (Gaussian) standard deviations for pul- 
sations at ~4.7 s could be obtained only by applying the FFT to 
the 1.9-4.9 keV flare events and using a time bin of Is. Different 
choices of the time binning would led to a result comparable to 
that obtained with our Fourier analysis (detection signicance 3 
sigma). 

Since pulsations in the ASCA data seems to be detected with 
a suciently high statistical significance only when a specie en- 
ergy range is selected and only performing the analysis with the 
FFT technique and an ad-hoc time bin of 1 s, we conclude that 
this detection could be due to a statistical fluctuation. 



4. Siv/ff data analysis 

Pulsations i n the X-ray emission of IGR J18410-0535 were also 
reported by ISidoli et al] (l2008l) . These authors used the three 
Swift IXWT observations ID. 00030988001, 00030988004, and 
00030988005 (see Table |4]l. All data were collected in photon 
counting mode (PC, time resolution 2.5 s), and summed up to 
extract a single source photon event list and increa se the statis- 
tics. B y using a folding technique on this event list, ISidoli et aH 
(I2008h found a peak in the periodogram at 4.70088±0.0004 s 
with a reduced of 6.3 (9 degree of freedom). 

4.1. Timing analysis 

We reanalyzed the Swift XRT data of the above Swift-XRT ob- 
servations to verify the detection of pulsations . We us ed the same 
Swift data analysis described in ISidoli et aP (l2008h . XRT data 
were processed with standard procedures (xrtpipeline v. 0. 1 1 .6). 
Filtering and screening criteria were applied by using ftools and 
only event grades 0-12 were considered. We extracted a source 
event list (194 events) by using a circular region of 29 arc- 
sec centred on the best known position of IGR J18410-0535 
(iHalpern et al. ''2004) and a background event list (1 106 events) 
by using an annular region centered on the same position and 
with inner radius of 245 arcsec and outer radius of 630 arcsec. 
The radii of the background extraction region were chosen in or- 
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4.698 4.7 4.702 4.704 

Periodogram (s) 




4.698 4.7 4.702 4.704 

Periodogram (s) 



Fig. 16. Periodogram realized with the SwiftXKT data by using 
the source (upper panel) and background (lower panel) event 
lists extracted by summing up data in the three observations 
ID. 00030988001, 00030988004, and 00030988005. 



Table 4. Swift observations log of IGR J18410-0535. 



Obs. ID 


Instr. 


Stai't time 


Stop time 


Exp. 






(UTC) 


(UTC) 


(ks) 


00030988001 


XRT/PC 


2007-10-26 00:08 


2007-10-26 06:45 


1.4 


00030988004 


XRT/PC 


2007-11-05 09:08 


2007-11-05 09:28 


1.2 


00030988005 


XRT/PC 


2007-11-09 16:11 


2007-11-09 16:32 


1.3 



Notes. Exp. indicates the efi"ective exposure time of the observation. 



der to avoid any contamination from the tail of the PSF0. The 
arrival times of the source and background events were all con- 
verted to the SSB with the task barycorr. We applied the folding 
technique adopted bv lSidoli et al.l (1200 8) to the source and back- 
ground event lists, to search for a period of ~4.7 s. In the peri- 
odogram obtained from the analysis of source events (Fig. [T6l 
upper panel) the peak at 4.70088 iY^^Jd.o.f .-6.5/9) is visible. 
However, we noted that the same peak is also present in the peri- 
odogram obtained from the background event list (Fig. [16] lower 
panel). In this case, we obtained =23. The difference in the 
value of the x^^^ (a factor ~4-5) is consistent with the different 
quality of the statistics content of the two data samp les. We thus 
argue that the detection of pulsations reported by ISidoli et aP 
(l2O08h was merely due to an instrumental effect present in the 
photon arrival times of the SwftpCRT data. Investigations of 



^ See also |http://heasarc.gsfc.nasa.gov/docs/swift/analysis/xrt-Sw| 
guide_vl_2.pdf. 



similar effects in other Swift PiKT data-sets are currently going 
on and will be reported in the relevant on-line pagefl 



5. Discussion 

During the 45 ks-long XMM-Newton observation presented 
here, IGR J18410-0535 underwent a bright X-ray flare lasting 
about 15 ks. The rise of the flare occurred 5 ks after the begin- 
ning of the observation, thus we were able to follow in detail the 
whole event to the return to quiescence. A time-resolved spec- 
tral analysis of the observation revealed that the source X-ray 
emission could be described well in the different selected time 
intervals by using an absorbed power-law model. The photon 
index of the power law remained virtually constant at least for 
the first 20 ks of the observation, but in this period the absorp- 
tion column density underwent dramatic changes. We measured 
a significant rise in the A^h from ~3 to 20x10^^ cm"^ across the 
transition from quiescence to the peak of the flare, followed by 
a sudden decrease to a value of 10x10^^ cm"^ which then re- 
mained stable for the following 4 ks. In the subsequent 1 1 ks, 
the Ah progressively increased to a value of ~50xlO^^ cm"^, 
and the spectral properties of the source changed dramatically 
during the last 23 ks of observation. In this time interval, we 
measured an abrupt flattening of the source power-law photon 
index and revealed the appearance of a prominent emission line 
at ~6.4 keV. Similar behavior is usually observed during the oc- 
currence of X-ray eclipses in some of th e sgHMXBs and in only 
one SFXT (IGR J16479-4514, see e.g.. lBozzo et"al]|2008bt and 
references therein). When the NSs in one of these systems is 
obscured by its supergiant companion, the X-ray emission pro- 
duced by the accretion process close to the surface of the com- 
pact object is progressively absorbed in the photosphere of the 
supergiant star Only the X-rays reflected by matter in the stellar 
wind remain visible to the observer and a fluorescence emission 
line at -6.4 keV becomes very prominent. 

In the present case, that IGR J18410-0535 underwent a 
bright X-ray flare after which the source turned off sharply be- 
cause of an X-ray eclipse appears to be unlikely. Below we dis- 
cuss an alternative scenario, in which the entire event observed 
by XMM-Newton is interpreted as the accretion of a massive 
clump onto the compact object hosted in IGR J18410-0535. 

5.1. Ingestion of a massive clump 

The idea that dense "clumps" of matter might be generated in 
the wind of supergiant stars as a consequence of different mag- 
netohydrodynamic instabilities has found some su pport in the 
results of several numerica l simulations (Runacr es & Owockil 
120021 lOskinova et al.ll2007l and referen ces therein) and in ob- 
servations of isolated OB stars (see e.g..'Eversberg et al]|1998l: 
Martins 2011). Even though the physical parameters of these 
clumps (i.e., mass, radius and density) are still largely unknown 
because of poorly estimated theoretical and observational pa- 
rameters, a clumpy wind scenario is the founding hypothesis 
of most models developed so far to interpret the X-ray variabil- 
ity of the SFXTs {i n't Zand 2005; Walter & Zurita Heras 2007|; 
IBozzo et"an l2008a[ and references therein). All these models 
predict that fast X-ray flares can be produced by the sporadic 
capture and accretion of one of these clumps onto the NSs hosted 
in the SFXTs. The role of these dense clumps distributed all 
around the NS orbit is two-fold. Clumps simply passing in front 



' http://www.swift.ac.uk/xrtdigest.shtml. 
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of the X-ray source cause dimming or even obscuration, and dis- 
play the signatures of photoelectric absorption, at least during 
the ingress and egress stages. In addition to these phenomena, 
clumps that lead to increased accretion also give rise to large 
variations in the X-ray luminosity. To reach the required X-ray 
luminosity level, different densities of the clumps might be re- 
quired, depending also on the strength of the NS magnetic field 
and the value of its spin period. Both the rotation of the NS and 
the extent of its magnetosphere can act as gates that halt most 
of the accretion during quiescence and release it during the out- 
burst, thus signi ficantly altering the total dynamic range of the X - 
ray luminosity dGrebenev & Sunvaevl2007l : lBozzo et al .I2008ah . 

In the case of IGR J18410-0535, the entire flaring episode 
observed by XMM-Newton might be due to the accretion of a 
single massive clump. 

- f=0-3200 s (time interval A): Pre-flare quiescence. At 
the beginning of the XMM-Newton observation, the 
source was caught at a relatively low X-ray flux 
(3.2x10"'"^ erg cm"^ s~'). For a source distance of 
3 kpc, this would translate into an X-ray luminosity of 
3.5x10-'^ erg/cm^/s, i.e. a value typical of SFXTs in qui- 
escence. The spectral parameters, A^h, T, that we measured 
in this time interval (spectrum A in Table [1]) are compati- 
ble with those inferre d for the "very low state" reported by 
iRomano et a D (l2009l) . The relatively low flux and absorp- 
tion column density measured in this state with respect to the 
other average X-ray emission states ("medium" and "high") 
suggests that the NS at the beginning of the XMM-Newton 
observation might have been located in a tenuous region of 
the stellar wind such that accretion was taking place only at 
a very low level or that it was inhibited by magnetic and/or 
centrifugal barriers. 

- f=3200-5220 s (time intervals B, C, and D): Rise of the flare. 
Around f =3200 s, the XMM-Newton lightcurve showed a re- 
markable rapid increase in the count-rate from IGR J 184 10- 
0535. This is firstly noticeable in the hard energy band (4- 
12 keV, see Fig.|2]), thus indicating a high level of absorp- 
tion. During the rise of the flare, we measured a significant 
increase in A^h, reaching a value of ~2xl0^^ cm"^ at the 
top of the flare (f^5200 s). At this stage, the source X-ray 
flux increased to 3.3x10"'" erg cm"^ s"', corresponding to 
an X-ray luminosity of 3.6x10-'^ erg/s. A similar simulta- 
neous increase in the X-ray luminosity and absorption col- 
umn density can be well understood by assuming that a mas- 
sive clump approached the NS magnetosphere and filled the 
immediate surroundings with high density material. The in- 
crease in the local density can easily lead to a significant 
increase in the A^h along our line of sight and eventually 
compress the NS magnetosphere to open the magnetic and/or 
centrifugal barrier and permit direct accretion onto the com- 
pact object. In this interpretation, we can use the observa- 
tional results in Table[T]to estimate the physical properties of 
the clump. 

In wind accretion theory, the typical time scale on which 
matter is accreted and reaches the surface of the NS from the 
magnetospheric boundary is of the orde r of the local free-fa ll 
time, i.e. hundreds of seconds (see e.g. jBozzo et al ] l2008ah . 
The duration of the flare observed from IGR J18410-0535 is 
~ 15 ks, thus this time is likely to be linked to the radial extent 
of the clump. Therefore, we can estimate 

T^d - l/2VwAfflare = 8 X 10"v„8 cm, (1) 

where Rd is the radius of the clump and v^s is the relative 
velocity between the clump and the NS in units of 10*^ cm/s 



(we neglect here the orbital velocity of the NS and consider 
a spherical clump moving with the same velocity as the sur- 
rounding stellar material; Lepine & Moffat 2008). To a first 
order approximation, only the material that falls inside the 
accretion radius, R^^c, of the NS is accreted, and thus 



Md = {R,ilR.,,,yM.,, 



(2) 



where R.^^^-2GMnsl\%,, Mns = 1-4Mq is the mass of the NS, 
and Mace is the total mass accreted during the flare. The latter 
can be estimated by integrating the X-ray flux measured dur- 
ing the flare as a function of time in Fig.|7]and using the rela- 
tions Funahs^LxliAncf) and Lx-GM^sM^cJRns (where d is 
the source distance, /?ns=10^ cm is the NS radius and M-^^c 
is the mass accretion rate). We found Macc=l-5xlO'^c/3kpc g 
(where dikpc is the source distance in units of 3 kpc) and thus 



Md = 1.4xl02V^,4p,g 



(3) 



which clearly has a strong dependence on the wind velocity. 
Another equation that can be used as a cross-check of these 
results is that relating the mass and radius of the clump to the 
expected absorption column density caused by its presence 
in the vicinity of the NS, i.e. 



A^H - Md/(/?>ip) = 1.3 X lO^^vts^/'kpcCm 



-2 



(4) 



As the A^H measured by the spectral analysis at the top of the 
flare is ~2xl0^^ cm"^, Eq.|4]implies a wind velocity slightly 
lower than 10^ cm/s. We note that the estimated values of M^ 
and Rc\ are in qualitative agreement with those expected ac- 
cording to the clumpy wind model developed by Ducci et al] 
(j2009t). 

f=5220-6320 s (time interval E): Top of the flare. The results 
of the spectral analysis obtained during this time interval can 
also be clearly explained within the scenario depicted above. 
The time-resolved spectral analysis carried out in Sect. 12.21 
showed that the initial rise in the A^h, interpreted above as 
being due to a massive clump approaching the NS, was fol- 
lowed by a sudden decrease in the absorption column density 
(by a factor of ~2) when the system reached the highest lu- 
minosity. According to the calculation of Krolik & Kallm^ 
(|1984|) . a similar drop in A^h is to be expected because of 
the heating and photoionization of the clump material by the 
higher X-ray flux. In particular, a major change in the opacity 
of the material with respect to the X-ray photons is expected 
when the so-called ionization parameter H^IO. The latter is 
defined as 



(5) 



where Lis=Lx/(lO^^ erg/s) is the luminosity coming from 
the X-ray source, r^R^i is the distance of the material from 
the source, and T -^IQ^ K is the temperature of the clump 
dPucci et alj l2009l) . We therefore expect H to be close to 
the critical value for a luminosity similar to that reached by 
IGRJ18410-0535 at the peak of the flare. During this time 
interval, we also detected a significant iron line with a cen- 
troid energy of 6.56;^q q5 keV. This is somewhat higher than 
the Fe Ka line of neutral iron (~6.4 keV), thus suggests that 
the X-ray flux produced by accretio n onto the NS p artly ion- 
ized the clump matter (according to lKallman et al.ll200 4. the 
centroid energy we measured at the top of the flare requires 
the presence of iron ions with ionization stages higher than 
FeXXI). 
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f=6320-9520 s (time intervals F, G, H, and I): Beginning 
of the flare decay. During this interval, the source X-ray 
flux was beginning to decrease. When it dropped below 
~6xl0"" erg cmT^ s"', a new rise in the A^h is measured 
(spectrum I). This is to be expected according to our inter- 
pretation where the X-ray flux decreases below the level re- 
quired for photoionising the clump material. A further pos- 
sible indication of the change in the ionization stage of the 
material around the NS can be deduced by the energy of the 
iron line detected during interval I. This is somehow lower 
than that measured at the top of the flare and close to the 
value expected in case of neutral iron (see Fig.fTTTl. 
An enigmatic finding that emerged from the analysis of the 
XMM-Newton data during this time interval is the drop in the 
source count-rate visible in Fig. [1] around f=8500-9500. As 
we discussed in Sect. 12.21 during this interval no particular 
spectral change was observed in the X-ray emission (we note 
only a marginal increase in the A^h with respect to the spectra 
extracted from the nearby time intervals and detected a sig- 
nificant iron line at ~6.4 keV, see Table [TJ. A closer inspec- 
tion of the XMM-Newton lightcurve of the source around 
?=8500-9500 revealed that the apparent drop in count-rate 
occurred in both the low (0.3-4 keV) and hard (4-12 keV) 
energy bands and that several small flares took place during 
this period (see Fig.[3]i. We suggest here that a qualitatively 
similar variations in the X-ray flux from the source might 
be expected because of some instabilities taking place close 
to the NS magnetosphere (e.g., when the decreasing accre- 
tion rate of material from the clump is close to the threshold 
for the opening/closing of the centrifugal and/or magnetic 
barrier IBozzo et al. 2008a). Alternatively, it can be related 
to some small structures in the stellar wind located close or 
within the clump itself. 

f=9520-42840 s (time intervals J, K, L, M, N, and O): Late 
stages of the flare decay. The spectral analysis revealed a 
progressive increase in the absorption column density, that 
reached a value of A^h-50x10^^ cm"^ in the time interval N. 
Even though an increase of the A^h is still expected at these 
times because of the lower ionizing X-ray flux, the value 
above is significantly higher than that measured when the 
source reached the peak of the X-ray flux. This result can 
be explained by assuming that, after the occurrence of the 
flare, the clump was located in front of the NS along our 
line of sight to the source. This would agree with the dra- 
matic change in the X-ray spectrum of the source observed 
at the end of the flare during interval O (f= 18870-42840 s). 
We suggested in Sect. 12.21 that a reasonable fit to this spec- 
trum could have been obtained with a model comprising 
two power-law components with different absorptions plus 
an iron line at 6.4 keV. This model is usually adopted to in- 
terpret the X-r ay emission of sg HMXBs in ec lipse (see e.g, 
Ivan der Meer e t al. 2005; Torreio n et alJl2010l for a review), 
but in the case of IGR J18410-0535, we assume that the ob- 
scuration of the NS was caused by a massive clump rather 
than by its supergiant companion. The more absorbed power 
law represents the X-ray emission caused by the accretion 
onto the NS that is strongly extinguished by the presence of 
the clump. The less absorbed power-law component is in- 
troduced to take into account the scattered emission in the 
stellar wind material spread all around the binary system, 
and might also comprise the X-ra y emission produce d by the 
shocks in the wind itself (see e.g.. Bozzo et al.ll2010l and ref- 
erences therein). This emission is unaffected by the presence 
of the clump and is thus only seen through an absorption col- 



umn density that is compatible with the Galactic value. The 
iron line is caused by X-ray irradiation of cold iron in the 
wind of the supergiant star (see, e.g. Kallman et al. 200%, 
and its relatively high EW also support the idea that most 
of the continuum emission is suppressed by the obscuration 
of the clump. However, if we consider that the obscuration 
of the X-ray source is taking place during the entire time 
interval O, some ad hoc assumptions on the system geome- 
try would probably be required to explain how a clump with 
the physical dimensions estimated earlier in this section and 
moving with a velocity of the order of ~10^ cm/s can hide 
the NS for about ~20 ks. We suggest instead that the obscu- 
ration occurred only during the first ~3 ks of the time inter- 
val O (i.e. time interval Oi). Once the clump moved away, 
the accretion onto the NS decreased substantially and the X- 
ray emission detected during the last ~20 ks of observation 
(time interval O2) was partly due to the residual accretion 
and partly to the scattering of these X-rays in the wind of 
the supergiant s tar and the shocks o ccurring within the wind 
itself (see e.g.. 'Bo zzo et aljF2010l and references therein). 
Even though the statistics of the data accumulated during the 
time intervals Oi and O2 was of too low quality to enable 
us to draw a firm conclusion, the analysis of these spectra 
carried out in Sect.|2]is consistent with this idea. We showed 
that a fit to the Oi spectrum required both the presence of 
the highly absorbed and the less absorbed power-law com- 
ponent, whereas the spectrum O2 could be described well 
using only the less absorbed power-law. However, given the 
very poor spectral information available during the time in- 
terval O2, we cannot presently rule out that after the clump 
moved away, the obscuration of the NS by its supergiant 
companion contributed to reduce the flux measured at this 
time (^9x10"'"^ erg cm"^ s"\ 1-10 keV, i.e a factor of ~3.5 
lower than that measured in the interval A, see Table [T]i. 
As a final remark, we also note that an iron line with a cen- 
troid energy and wEW compatible with that measured dur- 
ing the time interval I and E by XMM-Newton was detected 
from IGRJ18410-0535 during the ASCA observations car- 
ried out in 1994 and 1999 (see Sect. [TJ. On those occasions, 
the line was clearly visible in the spectra of the source ex- 
tracted in the flux range (O.l-l.O)xlO"" erg cm"^ s"' (2- 
10 keV), i.e. significantly different from those corresponding 
to the XMM-Newton detections (see Table[T]l. No convincing 
indication of the presence of an iron line could be found in 
the XMM-Newton spectra extracted at similar flux levels (in- 
tervals L, M, N; see also Sect. I2.2l i. As the centroid energy 
of the iron line, together with its EW and normalization, pro- 
vides information on the ionization state of stellar wind ma- 
terial and the amount of material around the NS, these results 
support the idea that the conditions of the stellar winds in the 
SFXT so urces can change significantl y during the NS orbits 
(see e.g.. lWalter & Zurita Herai2007h . 

5.2. IGR J1 841 0-0535: another SFXT hosting a 
slow-spinning NS? 

The timing analysis of the XMM-Newton data did not reveal any 
pulsations at 4.7 s, which is the periodicity detected during pre- 
viously ASCA and Swift observations of IGR J18410-0535 (see 
Sect. [TJ. The relatively tight upper limits on the pulsed fraction 
derived from the XMM-Newton data (see TableO, convinced us 
to reanalyze the archival ASCA and Swift data where pulsations 
were detected. In the first case, our results revealed that the de- 
tection might have been due to a statistical fluctuation, while in 
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the second case the detection was most likely related to an in- 
strumental problem. In addition the analysis of the latest Swift 
observations performed during the outburst of the source that 
occurred on 2010 June 5 could not confirm the presence of pul- 
sations ( Romano et al 12010b ). These results suggest that the NS 
in IGR J18410-0535 might not be pulsating at a period of 4.7 s. 

Among the ~15 sources in the SFXT class, so far only 
IGRJ11215-5952, IGR J 18483-03 11, and IGRJ16465- 
4507, displayed unambiguously cohere nt pulsations at 
-186.78 s, -21 s, and -228 s, respectively (ISidoli et alJl2007t 
IWalter & Zurita HerasI l2007t iGiunta et al.l l2009h . However. 
IGR Jl 12 15-5952 proved to be a very peculiar source, as it is the 
only SFXT displaying regularly periodic outbursts connected 
with the NS passage at the periastron of the system, while 
IGR J 18483-03 11 and IGR J16465-4507 were classified as "in- 
termediate SFXTs" due to their somewhat longer outbursts (up 
to a few days as opposed to a few hours) and sma ller dynamic 
range in the X-ray fl ux vari ation (<10-^, see e.g. Rahoui et aLj 
[2008; R omano et al] l2010cl: iLaParola et all 12010: Clark etalJ 
[201 0). All SFXTs displaying a very large dynamic range in the 
X-ray luminosity (>10'*-105), including IGR J18410-0535, seem 
instead to be "non-pulsating" sources (see the cases of, e.g. 
XTE J 1739-302 and IGRJ08408-450 3 and IGR J 17544-26 19 
iBozzo et"ar I I2OIOI: iRampv et al.ll2009l) . In some of these cases, 
pulsations were searched deeply up to periods of seve ral hun- 
dred s econds, and no evidence for them was found (Sm ith et alJ 
I2006t iBozzo etal J 120091 l2010h . These non-detections are stiU 
consistent with the possibility that these SFXT sources might 
host NSs with very long spin periods (>1000 s) and that their 
extreme X-ray variability co uld be rel ated to the centrifugal and 
magnetic gating mechanism dBozzo et aljr2008ah . 



6. Conclusion 

We have reported on the detection with XMM-Newton of a rela- 
tively bright flare from the SFXT source IGR J18410-0535. The 
lightcurve and spectral analysis presented in Sect. |2l provided 
several indications that the event was entirely due to the accre- 
tion of a massive clump onto the NS hosted in this system. Even 
thought the occurrence of an eclipse by the companion star to- 
ward the end of the observation could not be completely ruled 
out (see Sect. |5]l, we provided convincing evidence in favor of 
the "ingestion of a clump" scenario. We have shown that: 

- the rise of the flare was followed by a sudden increase in 
the absorption column density in the direction of the source. 
This is expected if the event is triggered by the presence of a 
massive clump approaching the NS. 

- at the peak of the flare, a significant drop in the absorption 
column density occurred; this was ascribed to the effect of 
the X-ray photoionization caused by the higher X-ray lumi- 
nosity of the source (this is also supported by the evidence of 
a change in energy of the iron line centroid with the source 
luminosity). 

- the gradual increase in the absorption column density in the 
direction of the source during the decay from the top of flare, 
together with the detection of a prominent iron line in the last 
part of the XMM-Newton observation, could be interpreted 
as implying that the remains of the clump first covered the 
X-ray source for a few ks (-3 ks) and then moved away 
from the NS. The residual X-ray emission observed in the 
last -17 ks of the observation might be due to some residual 
accretion taking place onto the NS. 



In accordance with this interpretation, we have provided a rough 
estimate of the mass and radius of the clump by assuming a 
spherical shape. We note that, because of our poor knowledge 
of the properties of the binary system IGR J18410-0535 and the 
lack of a clear confirmation of its spin period (see Sect. 15.21 ). it 
was not possible to take into account in these calculations the 
effect of the NS magnetic field according to the model devel- 
oped by Bozzo et al. (2008b). The lack of clear detection of rel- 
atively short spin periods in all the most extreme SFXT sources, 
including IGR J18410-0535, is still consistent with the possibil- 
ity that these object might host NSs with very long spin peri- 
ods (>1000 s) and that their extreme X-ray variability could be 
partly related to the centrifugal and magnetic gating mechanisms 
( iBozzo et al.ll2008ah . 
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Fig. 17. An artist's view of the "ingestion of a clump" described in Sect. 15.11 Different elements in the figure are not represented to 
scale. The figures illustrate qualitatively the main phases of the event in chronological order from the top left to the bottom right. In 
each figure, we plot in the top right corner the XMM-Newton lightcurve of the event already shown in Fig. [7] and highlight with a 
red color the data points in the plot that correspond to the physical scenario described in the picture. From top left to bottom right: 
(a) at the beginning of the XMM-Newton observation the NS is accreting from a low density material, such that the observed X-ray 
luminosity is relatively low (~4xl0^^ erg,?"'); (b) the NS encounters the clump and XMM-Newton observes a rapid rise in the X-ray 
flux from the source; (c) the higher X-ray luminosity (-4x10^^ ergi"') causes the photoionization of the surrounding material and 
a decrease in A^h; (d) accretion decreases and the source undergoes a rapid decay in the X-ray flux; (e) the remains of the clump 
move in front of the NS along the observer line of sight and obscure the X-ray source; (f) the remains of the clump move away and 
the source returns to its quiescent state (~10^^ ergs"', some residual accretion might still take place). 
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